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Abstract 

A  recently  developed  ternary  phase- field  model  [Wen  et  al.,  Acta  Mater,  submit¬ 
ted]  is  applied  to  study  the  precipitation  behavior  in  Ni-Al-Cr  alloys.  The  effect 
of  diffusivity  and  initial  compositional/ordering  profiles  on  the  growth  process  of 
precipitates  are  examined.  Specifically,  the  composition  evolution  with  time  in  the 
precipitates  (7')  is  simulated  and  compared  with  experimental  observations  by  Sud- 
brack  et  al.  [Acta  mater.,  54:3199,2006].  The  simulation  work  reveals  that  the  Cr- 
trapping  phenomenon  can  only  be  reproduced  with  no  significant  partitioning  of  Cr 
at  the  nucleation  stage.  By  looking  into  the  composition  dependent  diffusivity  and 
thermodynamic  driving  force  for  both  Al  and  Cr,  it  is  found  that  the  key  differ¬ 
ence  resides  in  the  much  stronger  composition  dependence  of  the  diffusivity  of  Al. 
The  diffuvisity  of  Al  is  enhanced  by  roughly  an  order  of  magnitude  across  regions 
near  or  at  a  7/7'  interface  region.  On  the  contrary,  the  diffusivity  of  Cr  does  not 
show  a  strong  composition  dependence.  This  is  believed  to  be  the  major  factor  that 
has  caused  Cr  being  left  ’trapped’  during  the  growth  of  the  precipitate  while  Al 
approaching  to  the  equilibrium  concentration  quickly  to  minimize  the  total  energy. 

Key  words:  ageing,  Ni-Al-Cr  alloys,  ternary  phase-field  models,  precipitation, 
Cr-trapping 


1  Introduction 


The  compositional  evolution  inside  the  7r-precipitate  and  7-matrix  was  re¬ 
cently  measured  with  atom-probe  tomography  for  Ni-5.2  Al-14.2  Cr  at.%  spec¬ 
imens  aged  at  873K  [1],  The  results  were  analysed  by  comparing  with  CAL- 
PHAD  assessments.  The  measured  values  were  largely  in  line  with  CALPHAD 
predictions.  For  example,  the  experimental  values  for  the  equilibrium  com¬ 
position  of  the  7-matrix  are  near  the  CALPHAD  predictions.  The  observed 
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decreasing  concentration  of  A1  in  the  interior  of  the  y'-precipitate  with  time 
followed  the  same  trend  as  predicted  by  CALPHAD.  The  most  striking  dis¬ 
crepancy  was,  however,  the  evolution  of  Cr  concentration  in  the  y'-precipitate. 
The  observed  concentration  of  Cr  in  the  precipitate  decreases  monotonically 
with  the  growth  of  the  precipitate  while  the  CALPHAD  approach  predicts  the 
opposite.  This  phenomenon  was  referred  to  as  Cr-trapping.  The  objective  of 
this  work  is  to  understand  the  mechanism  behind  the  Cr-trapping  phenomenon 
using  our  recently  developed  ternary  phase-field  model  [2], 

The  observed  trend  in  Cr  concentration  in  the  precipitate  was  attributed  to 
the  kinetically  trapped  Cr  atoms  [3]  based  on  a  mean  field  evaluation  that  the 
intrinsic  diffusivity  of  Cr  is  slightly  smaller  than  that  of  A1  [4],  This  difference 
in  diffusivity  is,  however,  not  significant  enough  to  be  the  sole  factor  to  cause 
such  a  dramatic  change  in  the  trend  of  the  Cr  evolution  in  the  y'-precipitat. 
Sudbrack  et  al  [3]  reported  the  presence  of  Ni3Al-type  Ll2  short  range  order¬ 
ing  domains  in  the  as-quenched  alloy,  which  was  most  likely  resulted  from  fast 
Al  diffusion  during  the  quench  from  the  solutionizing  temperature  of  1123K. 
It  is  speculated  that  the  initial  inhomogeneity  in  structure/composition  may 
have  played  an  important  role  in  the  nucleation  of  y'-precipitate  resulting  in  a 
change  in  the  evolution  behavior  of  Cr  with  the  growth  of  the  precipitate.  An¬ 
other  possible  factor  is  the  thermodynamic  driving  force  considering  that  the 
diffusion  is  controlled  by  the  product  of  the  diffusion  mobility  and  the  gradient 
of  the  thermodynamic  driving  force.  The  driving  force  might  be  much  smaller 
for  Cr  diffusion  as  compared  to  that  for  Al  diffusion.  As  such,  the  Cr-trapping 
could  be  mainly  a  thermodynamic  effect  rather  than  a  kinetic  one.  In  summary, 
the  Cr-trapping  could  at  least  be  attributed  to  three  factors,  i.e.  smaller  dif¬ 
fusivity  of  Cr,  some  kind  of  favorable  initial  nucleus  structural/compositional 
profiles,  and  smaller  thermodynamic  driving  force  for  Cr  diffusion.  It  should  be 
pointed  out  that  these  three  factors  are  interrelated.  For  example,  the  initial 
nucleus  structural/compositional  profiles  may  be  largely  determined  by  defect 
structure  of  the  solute  solution,  which  is  strongly  affected  by  the  specific  heat 
treating  of  the  sample,  e.g.  the  cooling  rate.  But  it  can  also  be  controlled  by 
the  thermodynamic  driving  force  and  diffusivity  of  Al  and  Cr.  The  focus  of  this 
work  is  to  identify  the  most  dominant  factor  that  has  caused  the  Cr-trapping 
phenomenon.  In  what  follows,  the  ternary  phase-field  model  is  described  and 
some  simulation  results  and  a  brief  discussion  are  presented  in  Sec. 3. 


2  The  ternary  phase-field  model 


A  more  detailed  description  of  our  ternary  phase-field  model  was  given  else¬ 
where  [2],  The  model  linked  directly  to  commercial  CALPHAD  software.  A 
diffusion  mobility  database  was  implemented  that  included  the  contribution 
from  the  ordered  state  allowing  us  to  have  a  better  description  of  diffusion 
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behavior  near  the  7/7'  interfaces  and  inside  the  ordered  7'  phase.  The  mis¬ 
fit  effects  of  the  partitioning  of  the  two  solute  elements  were  also  taken  into 
account. 

Six  long-range  order  ( Iro )  parameters,  i.e.  rfc  (r,  t)  and  two  concentration  field 
variables  xn  (r,  t)  are  used  to  describe  the  two-phase  microstructure  (7  +  7') 
in  the  Ni-Al-Cr  systems  [2],  Here  k  =  1,2,3  are  the  three  components  of 
Iro  parameters,  n  =  1,2  representing  A1  and  Cr,  respectively,  r  and  t  are 
spatial  coordinates  and  time.  These  eight  field  variables  can  fully  describe  the 
concentration  inhomogeneity  of  solute  A1  and  Cr  in  the  two  phases  and  their 
site  fractions  in  the  ordered  domain  structures  in  7'. 

The  microstructural  evolution  can  be  simulated  by  solving  the  time-dependent 
Ginzburg-Landau  equation  for  the  six  Iro  parameters  and  using  the  non-linear 
Cahn-Hilliard  diffusion  equation  for  the  two  concentration  fields: 

(1) 
(2) 

where  L  is  the  structural  relaxation  constant,  Mij  are  the  diffusion  mobilities, 
Vm  is  the  molar  volume,  and  F  is  the  total  free  energy  of  the  system. 

For  the  coherent  precipitation  under  consideration,  the  total  free  energy  con¬ 
sists  of  the  elastic  strain  energy  (Fei)  and  the  chemical  free  energy  ( Fch ): 


drjj 

dt 

Ox, 


=  -L 


5F 


dt 


m  ' 


F  —  Fei  +  Fch  (3) 

The  elastic  strain  energy  arises  because  of  the  lattice  mismatch  between  the  7 
and  7'  phases.  The  misfit  strain  is  assumed  to  depend  on  the  compositions  of 
the  two  solute  elements  in  the  two  phases.  The  elastic  energy  of  a  two-phase 
mixture  can  then  be  calculated  as  a  function  of  the  concentration  profiles  [2]: 

V"  _  2  _ 

Fei  G Cijkl^ij  2  j  £ kl  ip')  t-t^Xp  (4) 

1  p=  1 

v  2  2  _ 

+  E  E  e°ij  (p)  eki  («)  AxP' AxQ 

Z  p= 1  <7=1 

1  _  ^  2  nt  (]3g- 

-oEE  /  7^3 BPQ  (n)  {A^}g{Ax  Jg, 

Zp=lq=lJ  (27T) 

where  ^4^)  is  the  average  volume  of  (A),  V  is  the  total  volume  of  the  system, 
Cijki  is  the  elastic  moduli  tensor, (p)  is  the  stress  free  transformation  strain 
tensor  associated  with  the  partitioning  of  solute  p,  Axp  is  the  deviation  of  the 
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local  concentration  of  solute  p  from  its  bulk  concentration.  The  notation  J' 
indicates  that  the  point  n  =  0  is  excluded  from  the  integration,  and  Bpq  (n)  is 
a  two-body  interaction  potential  given  by: 

Bpq{n )  =  ma-j  ( p )  Qjk  (n)  okl(q)nh  (5) 

where  n  =  jJj  is  a  unit  vector  in  reciprocal  space  and  n*  is  the  ith  component, 
a°j  =  Cijki^h  ( p )  and  (n)  is  a  Green  function  tensor  which  is  inverse  to 
the  tensor  fi”1  (n)  =  i,  {Axp}s  is  the  Fourier  transform  of  Axp,  and 

{A xPYg  is  the  complex  conjugate  of  {Axp}  . 

The  chemical  free  energy  consists  of  the  local  chemical  free  energy  and  the 
gradient  terms  arising  from  compositional  and  structural  inhomogeneities  [5], 
i.e., 


c  1 

Fch  =  I  f(xAhXCr,rfAi,rfCr,T)  + -  ^  (n^VxiVxA 

JV  ^  A  A—  A1  /"tv. 


i,j-Al,Cr 


+\Kn  E  E(Vr^)5 

Z  i=Al,Cr  P 


dV 


(6) 


where  f  (xAi,  %cr ,  rf/u,  Vcn  F)  is  the  local  chemical  free  energy  that  defines  the 
basic  thermodynamic  properties  of  the  system,  k,x.  .  and  Kq  are  gradient  energy 
coefficients,  and  it  is  understood  that  p  =  1,2,3  for  the  three  components  of 
the  Iro  parameters. 


In  the  CALPHAD  method,  the  local  chemical  free  energy  is  expressed  as  a 
sum  of  a  reference  free  energy,  an  ideal  mixing  free  energy,  and  excess  mixing 
free  energy: 


G  (r,  T)  =  Gref  (r,  T)  +  Gld  (r,  T)  +  Gex  (r,  T) .  (7) 

The  three  energy  terms  depend  on  the  site  fractions  of  the  two  solute  elements 
as  follows: 


cref  (r^)  =  EEEE  vl  (r)  y)  (r)  vl  (r)  yf  (r)  Gm  CO ; 

i  1  j  1  k  1  fc=a 

(r,  T)  =  RT  ±  f  £  y‘H  (r)  In  [jJ  (r)] ;  (8) 

s—  1  n— 1 

4  4  TV  TV  Nipq 

G’*(T,T)  =  '£'£'£ y’p W E < W KM- E K M - < wrc?, (m, r> , 

s=l  t^sp=l  g=l  m= 1 

where  N  is  the  total  number  of  elements  in  the  alloy,  and  yf  (r)  is  the  site 
fraction  of  element  i  on  sublattice  s  with  normalizations  Vp  (r)  =  1  and 
4'Z)s=i  Up  (r)  =xp,  that  is,  summing  the  fractional  occupancies  over  all  the 
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alloying  elements  on  each  sublattice  must  equal  one,  and  summing  the  frac¬ 
tional  occupancies  over  all  sublattices  for  each  alloying  element  must  equal 
the  mole  fraction  of  that  element.  The  functions  G°]kl  (T)  and  G**  (m,  T)  are 
empirically  fitted  parameters  tabulated  as  functions  of  temperature  in  thermo¬ 
dynamic  databases.  R  is  the  gas  constant,  fs  is  the  site  fraction  of  sublattice 
s. 


The  local  chemical  free  energy  in  Eq.  7  is  expressed  as  a  function  of  the 
site  fractions  of  the  elements  on  each  of  the  four  sublattices.  However,  in  the 
phase-field  description  the  local  chemical  free  energy  is  a  function  of  local 
concentrations  and  Iro  parameters.  The  two  set  of  parameters  are  related  to 
each  other  as  described  below: 


Vn  Ob  t) 

vl  Ob  t) 

Vn  Ob  *) 
Vn  Ob 


xn  (r,  t) 

1  +  vl  Ob  t)  +  vl  Ob  t)  +  vl  (d  t) 

xn  (r,  t) 

1  -  vl  Ob  t)  -  vl  Ob  t)  +  vl  (d  t) 

xn  (r,  t) 

1  -  vl  Ob  t)  +  vl  Ob  t)  -  vl  (r,  t) 

xn  (r,  t) 

[1  +  vl  (r,  t)  -  vl  (r,  t)  -  vl  (r,  t)] 

(9) 


We  have  chosen  to  link  our  phase-field  code  to  the  commercial  CALPHAD 
software  Pandat  (CompuTherm,  LLC).  The  coupling  between  the  phase-field 
simulation  and  the  commercial  thermodynamic  database  involves  the  evalua¬ 
tion  of  the  partial  derivatives  of  the  local  chemical  free  energy  with  respect  to 
the  eight  site  fractions  (y*  (r,  t))  by  Pandat.  These  partial  derivatives  are  then 
used  to  evaluate  the  partial  derivatives  of  the  local  chemical  free  energy  with 
respect  to  the  phase-field  parameters,  i.e.  concentration  parameters  (xn  (r ,t)) 
and  the  Iro  parameters  (vl  (r,  t) ,  vl  (r,  t) ,  vl  (r,  t )),  through  the  chain  rules: 


df_ 

dxp 

df_ 

dr]$ 


Y  y  dG  dyn  1  23. 

s=i 

4  2 


zi  .1 


(10) 


The  diffusion  mobility  MtJ  varies  with  composition  and  structure.  Due  to  very 
limited  data  on  the  diffusion  mobility  of  the  ordered  phase,  most  simulations  to 
date  have  used  the  diffusion  mobility  of  the  disordered  fee  phase  to  represent 
that  of  the  overall  behavior  of  the  alloys.  In  other  words,  the  contribution 
from  the  ordered  phase  is  totally  ignored  although  it  is  well  understood  that 
the  diffusion  behavior  and  therefore  the  closely  related  diffusion  mobility  in 
the  ordered  phase  are  quite  different  from  that  in  the  disordered  fee  phase. 
Most  recently,  Campbell  has  developed  a  diffusion  mobility  database  that  is 
capable  of  describing  the  diffusion  mobilities  in  the  7'  phase  in  the  Ni-Al-Cr 
system  [6].  In  this  work,  we  adopted  this  database  for  the  description  of  the 
diffusion  mobilities  in  the  Ni-Al-Cr  system. 
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The  diffusion  mobility  is  expressed  as  [7] 

Mmai  =  T^XAI  [xcrXAlBcr  +  (1  -  xAif  BAl  +  xMxNiBNi  1  ;  (11) 

MAlCr  =  T^XaiXct  [~  (1  “  X Al)  BAl  ~  (1  —  XCr )  BCr  +  XNiBNi]  ; 

McrCr  =  T^Xcr  \xAlXCrBAl  +  (1  ~  Xcr)2  Bcr  +  , 

*m 

where  Bi  is  the  atomic  mobility  of  species  i.  The  atomic  mobility  is  usually 
expressed  as: 

B,  =  ^exp(-I\Q'/RT),  (12) 

where  R  is  the  gas  constant  and  Q*  is  the  total  activation  energy  of  species  i 
that  consists  of  contributions  from  the  disordered  (A  Qfs)  and  ordered  (A  Q°rd) 
states  as  described  below, 


A  Q*  =  A  Qfs  +  A  Q°rd.  (13) 

The  contribution  from  the  disordered  state,  A Qfs,  and  its  composition  de¬ 
pendence  can  be  obtained  based  on  Engstrom  and  Agren’s  work  [8].  The  con¬ 
tribution  from  the  ordered  state  is  described  below  as  in  Campbell’s  work  [6], 

=££  +  £  £  £  A  Q’SMvfli  -  w*] 

i  j^6i  i  j  k 

+  E  E  E  AQk:ij  [Vi  Vj  Vk  -  XiXjXk\ ,  (14) 

i  j  k 

where  the  index  m,  i,  j,  k  =  1,  2,  3  representing  the  three  components  Al,  Cr,  Ni 
in  the  ternay  system,  yf  is  the  site  fraction  of  component  i  on  the  a  sublattice 
that  can  be  readily  converted  from  the  current  4-sublattice  formulation  [9]. 
A Q™j,  AQ™.k,  A Q™{j  are  the  contributions  to  the  activation  energy  for  compo¬ 
nent  m  as  a  result  of  the  chemical  ordering  of  the  i,  j,  and  k  atoms  on  the  two 
sublattices.  These  parameters  are  well  documented  in  Campbell’s  work  [6]. 


3  Simulation  Results  and  Discussion 


As  mentioned  in  the  introduciton,  the  Cr-trapping  may  be  attributed  to  three 
factors,  i.e.  (a)  smaller  diffusivity  of  Cr,  (b)  some  kind  of  favorable  initial  nu¬ 
cleus  structural/compositional  profiles,  and  (c)  smaller  thermodynamic  driv¬ 
ing  force  for  Cr  diffusion.  Simulations  are  set  up  to  examine  how  these  three 
factors  affect  the  evolution  of  Cr  concentration  in  the  y'-precipitates  with  their 
growth.  Varying  diffusivity  from  benchmark  NIST  database  values  is  applied 
for  examining  factor  (a).  The  results  demonstrate  that  diffusivity  alone  can 
not  explain  the  Cr-trapping  phenomenon.  For  factor  (b),  some  representative 
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structural/compositional  profiles  need  to  be  chosen  to  study  their  effect  and 
some  rationale  is  given  below  for  the  choices.  The  contribution  of  the  thermo¬ 
dynamic  driving  force  is  analysed  in  terms  of  its  gradient  and  compared  with 
that  of  the  diffusivity  for  examining  factor  (c). 

For  a  binary  Ni-Al  system,  homogeneous  nucleation  is  considered  the  domi¬ 
nant  mechanism  for  7'-precipitate  nucleation.  Thermal  fluctuation  is  usually 
involved  in  the  formation  of  an  enriched  A1  embryo  with  LI  2-type  ordered 
structure,  which  serves  as  a  potential  nucleus.  In  numerical  simulations,  such 
a  nucleus  is  usually  represented  by  an  equilibrium  y'-precipitate  that  is  fully 
ordered  and  at  equilibrium  composition.  The  situation  becomes  more  complex 
in  a  multicomponent  system.  Take  the  current  ternary  system  Ni-Al-Cr  as  an 
example,  it  is  not  clear  whether  it  is  necessary  to  have  fully  ordered  A1  and 
Cr  simultaneously  and  at  their  respective  equilibrium  compositions  in  the  nu¬ 
cleus  for  it  to  survive.  Another  question  is  how  the  structural/compositional 
environment  surrounding  a  nucleus  (in  the  matrix)  such  as  a  partially  or¬ 
dered  matrix  will  affect  the  survivability  of  the  nucleus  and  its  growth.  To 
answer  these  questions,  it  is  necessary  to  study  the  effect  of  non-equilibrium 
structural/compositional  profiles  in  the  nucleus  and  matrix  on  the  growth  of 
precipitates  and  associated  Cr  concentration  evolution. 

In  our  phase-field  representation  of  the  two  phase  (  7  and  7')  microstructure, 
two  concentration  variables  and  two  sets  of  ordering  parameters  (for  A1  and  Cr 
respectively)  are  employed.  There  are  many  possible  combinations  as  initial 
profiles  for  the  precipitate  (7')  and  matrix  (7).  For  the  precipitate,  it  can  be 
assumed  with  different  degree  of  ordering  of  A1  and  Cr  for  testing  the  effect 
of  partial  ordering  or  different  degree  of  ordering  in  the  two  species;  and  var¬ 
ious  concentrations  of  A1  and  Cr  to  take  into  account  the  kinetic  effect  (e.g. 
smaller  diffusivity  of  Cr)  on  the  nucleus  formation.  While  the  matrix  can  rea¬ 
sonably  be  assumed  with  bulk  concentrations,  it  may  still  have  various  degree 
of  ordering  especially  for  A1  based  on  experimental  evidence  [3],  which  may 
affect  the  growth  kinetics  of  the  precipiate  and  associated  Cr  concentration 
evolution  with  its  growth.  In  this  work,  the  case  with  fully  ordered  A1  and 
Cr  and  equilibrium  concentrations  in  the  y'-precipitate,  along  with  disordered 
A1  and  Cr  and  bulk  concentrations  in  the  7-matrix  is  chosen  as  a  benchmark 
case.  This  case  is  referred  to  as  case  (i)  hereafter  for  convenience  of  further 
discussion.  Other  cases  examined  in  this  work  include  the  followings:  (ii)  fully 
ordered  A1  and  Cr  and  measured  7,-nucleus  composition  by  Sudbrack  et  al. 
for  the  precipitate  [1],  along  with  disordered  Al  and  Cr  and  bulk  concentra¬ 
tions  in  the  7-matrix;  (iii)  partially  ordered  (50%)  Al  and  Cr  and  equilibrium 
concentrations  in  the  7'-precipitate,  along  with  disordered  Al  and  Cr  and 
bulk  concentrations  in  the  7-matrix;  (iv)  fully  ordered  Al  and  Cr  and  equi¬ 
librium  concentrations  in  the  7'-precipitate,  partially  ordered  Al  (50%)  and 
disordered  Cr  and  bulk  concentrations  in  the  7-matrix;  (v)  homogeneous  and 
disordered  Cr  in  the  whole  simulation  cell,  fully  ordered  Al  and  equilibrium 
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concentrations  in  the  y'-precipitate  and  disordered  and  bulk  concentration  in 
the  7-matrix.  The  concentration  and  the  Iro  profiles  of  these  five  cases  are 
schematically  illustrated  in  Fig.  1  with  a  y'-precipitate  located  in  the  center. 

The  simulations  have  revealed  that  the  Cr-trapping  can  only  be  reproduced 
with  no  significant  partitioning  of  Cr  at  nucleation  stage.  This  is  a  major 
factor  that  contributes  to  the  Cr-trapping.  Another  factor  is  related  to  the 
composition-dependent  diffusivity  of  A1  and  Cr.  The  key  difference  resides 
in  the  much  stronger  composition  dependence  of  the  diffusivity  of  Al.  The 
diffuvisity  of  Al  is  enhanced  by  approximately  an  order  of  magnitude  across 
regions  that  have  significant  concentration  inhomogeneity  such  as  near  or  at  a 
7/7'  interface,  a  factor  that  is  ignored  by  mean  field  approaches.  In  contrast, 
the  diffusivity  of  Cr  is  not  at  all  sensitive  to  the  variation  in  concentration 
inhomogeneity.  This  is  another  major  factor  that  has  caused  the  Cr-trapping 
phenomenon. 

Fig.  2  shows  the  growth  of  a  single  precipitate  (top  row)  and  evolution  of 
Cr  concentration  with  time  in  the  precipitate  for  five  different  initial  struc¬ 
tural/compositional  profiles  described  above.  The  experimentally  measured 
Cr  concentration  evolution  is  also  included  for  comparison.  As  can  be  seen, 
the  Cr  concentration  in  the  first  four  cases,  i.e.  cases  (i)  through  (iv),  tends  to 
increase  with  time  with  growth  of  y'-precipitate,  which  is  the  opposite  of  the 
experimentally  observed  trend  as  shown  in  Fig.  2(d).  A  decreasing  Cr  concen¬ 
tration  with  time  is  observed  in  case  (v)  where  there  is  no  partitioning  of  Cr 
in  the  nucleus.  Similar  trend  is  also  found  with  the  growth  of  multiple  precipi¬ 
tates  (Fig.  3)  where  cases  (i)  and  (v)  are  included.  For  the  multiple  precipitate 
case,  nuclei  are  introduced  into  individual  cells  randomly,  but  with  the  prob¬ 
ability  related  to  the  corresponding  nucleation  rate  for  those  cells  [10]  [1 1]  [12] . 
These  results  suggest  that  the  experimentally  observed  trend  of  decreasing  Cr 
in  the  precipitate  with  its  growth  is  related  to  the  lack  of  partitioning  of  Cr 
in  the  nucleation  stage. 

It  was  pointed  out  that  a  slightly  higher  diffusivity  of  Al  relative  to  Cr  is  not 
significant  enough  to  be  the  sole  factor  to  result  in  the  decreasing  Cr  concen¬ 
tration  in  the  precipitate  with  its  growth.  To  validate  this  assessment,  some 
simulations  are  carried  out  by  varying  the  chemical  mobility  of  Cr  (Merer) 
(decreasing  for  the  purpose  of  exaggerating  the  difference  in  diffusivity  be¬ 
tween  Al  and  Cr)  that  in  turn  varied  (decreased)  the  diffusivity  of  Cr  since 
the  diffusity  is  proportional  to  the  mobility,  while  keeping  the  original  chem¬ 
ical  mobilities  of  Al  (m^ai)  aRd  the  cross  term  [Maict)  unaltered.  Fig.  4 
shows  the  simulated  evolution  of  Cr  with  the  chemical  mobility  of  Cr  set  at 
a  half  of  the  literature  value  [4],  which  means  that  the  difference  of  intrinsic 
diffusivities  of  Al  and  Cr  are  now  doubled.  This  may  account  for  the  possible 
disparity  between  experimental  value  and  thermodynamically  assessed  value 
due  to  gradient  effect.  As  can  be  seen,  the  mobility  (diffusivity)  change  slightly 
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affected  the  rate  of  evolution  of  Cr  for  the  five  cases  with  different  initial  nu¬ 
cleus  profiles  but  did  not  alter  the  trend  even  with  the  current  exaggerated 
difference  in  mobility  (diffusivity) .  This  proves  that  the  diffusivity  alone  can 
not  explain  the  Cr-trapping  phenomenon. 

To  gain  insight  about  the  respective  contribution  of  chemical  mobility  and 
thermodynamic  driving  force  to  the  Cr-trapping,  examination  of  the  two  key 
mobility  terms  {m^au  M^Cr^  and  the  gradient  of  the  thermodynamic  driving 

force  are  carried  out.  Fig.  5  shows  the  mobility  profiles  (top  row)  and 

the  ratio  of  M^Cr/MAl] Al  (bottom  row)  along  the  middle  cross  section  of  the 
simulation  cell.  Only  cases  (i)  and  (v)  (refer  to  Fig.  1)  are  included  in  the  chart. 
The  three  other  cases  present  almost  identical  profiles  as  the  case  started  with 
equilibrium  nucleus,  i.e.  case  (i).  The  results  show  a  significant  increase  in  A1 
mobility  in  the  core  of  the  precipitate  and  interface  region.  While  the  value 
of  MAllAl  is  slightly  higher  than  that  of  M^Cr  in  the  matrix,  the  difference  is 
greatly  increased  in  the  interface  region  especially  at  earlier  stages(Fig.  5(a)). 
In  case  (i)  with  equilibrium  nucleus,  the  A1  mobility  MAlAl  in  the  core  of 
the  precipitate  is  tripled  while  the  Cr  mobility  is  slightly  decreased.  In  case 
(v)  with  no  partitioning  of  Cr,  although  the  the  Cr  mobility  is  doubled  in 
the  core  of  the  precipitate,  the  A1  mobility  MAllAl  is  in  contrast  increased  by 
roughly  an  order  of  magnitude.  As  a  result,  the  ratio  of  the  two  mobilities 
{Mc*Cr/MAlAl)  in  both  cases  is  decreased  significantly  leading  to  a  relatively 
much  higher  diffusion  rate  of  A1  in  the  interface  region  (see  Fig.  5(c)).  With 
the  growth  of  the  precipitate  particle,  the  compositions  at  both  sides  of  the 
interface  approach  to  equilibrium.  As  a  result,  the  mobility  profiles  are  similar 
for  the  two  cases  at  later  stages  (Fig.  5(b)).  However,  the  enhancement  of  A1 
diffusion  in  the  interface  region  still  exists  as  evidenced  by  the  the  ratio  of  the 
two  mobilities  M^Cr/MAllAl  in  Fig.  5(d). 

Fig.  6  shows  the  gradient  of  the  thermodynamic  driving  force.  It  is  evident  that 
only  the  interface  region  has  shown  significant  gradient  of  the  thermodynamic 
driving  force.  The  magnitude  of  the  gradient  for  A1  evolution  in  the  interface 
region  is  comparable  with  that  for  Cr  evolution  in  the  whole  growth  process 
of  the  precipitate  as  represented  by  the  two  different  stages.  As  expected, 
the  gradient  of  the  thermodynamic  driving  force  near  the  interface  region  is 
greater  for  both  solutes  at  earlier  stages  and  it  decreases  monotonically  with 
the  growth  of  the  precipitate  and  the  system  is  approaching  to  equilibrium 
state  as  shown  in  Fig.  6(b). 

The  diffusion  rates,  i.e.  dxAi/dt  and  dxcr/dt ),  are  proportional  to  a  combined 
term  of  the  chemical  mobilities  and  gradient  of  thermodynamic  driving  forces 
as  shown  on  the  right  hand  side  of  Eq.  2.  The  diffusion  rates  are  shown  in  Fig.  7, 
which  represent  effective  driving  force  for  diffusions.  The  effective  driving  force 
exists  only  near  the  interface  region  where  a  mixture  of  two  phases  present.  It 
is  higher  in  value  in  the  earlier  stage  due  to  higher  degree  of  non-equilibrium 
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status  and  it  decreases  and  eventually  diminishes  at  later  stages  as  the  system 
approaches  to  equilibrium  status.  This  is  evidenced  by  the  two  snapshots  at 
two  representative  stages  in  Fig.  7(a)  and  (b).  Comparing  the  effective  driving 
force  for  A1  and  Cr,  it  is  significantly  higher  for  Al.  For  example,  the  effective 
driving  force  for  Al  is  approximately  one  order  of  magnitude  higher  than  that 
of  Cr  for  case  (v).  For  case  (i),  the  effective  driving  force  for  Cr  diffusion  is 
essentially  zero  while  that  for  Al  is  moderate.  It  should  be  pointed  out  that 
the  difference  in  effective  driving  force  for  Al  and  Cr  is  largely  due  to  the 
greatly  enhanced  diffusivity  in  the  interface  region  for  Al.  This  in  turn  leads 
to  a  much  faster  diffusion  of  Al  relative  to  Cr  and  leaving  the  Cr  behind  the 
interface  to  ’catch  up’  later,  thus  causing  the  Cr-trapping  phenomenon. 

Considering  that  the  much  enhanced  chemical  mobility  difference  in  the  7/7' 
interface  region  is  the  major  factor  for  the  Cr-trapping,  it  is  anticipated  that 
this  Cr-trapping  phenomenon  decreases  or  even  disappears  if  the  chemical 
mobility  difference  is  eliminated.  A  simulation  is  run  to  test  this  with  the 
chemical  mobility  of  Cr  set  equal  to  that  of  Al,  i.e.  M^Cr  =  MAl\t,  while  the 
cross  term  MA}Cr  keeps  unchanged.  This  has  in  effect  increased  the  chemical 
mobility  (diffusivity)  of  Cr.  The  results  are  shown  in  Fig.  8,  where  cases  (i)  and 
case  (v)  are  included.  For  comparison,  the  two  cases  shown  in  Fig.  2(c)  are  also 
included  as  references,  where  diffusion  mobilities  from  literature  databases  [6] 
are  used.  For  case  (i)  where  the  nucleus  has  equilibrium  composition  and 
equilibrum  Iro  configurations,  the  evolution  of  Cr  is  very  similar  to  that  in 
the  reference  case.  For  case  (v)  where  homogeneous  Cr  is  assumed,  a  rapid 
decrease  of  Cr  is  observed  in  the  first  few  hours  and  stays  flat  from  there.  In 
the  reference  case,  however,  a  steady  decrease  of  Cr  concentraion  is  observed  in 
the  whole  time  window  (approximately  300  hours  of  aging)  of  the  simulations. 
This  result  indicates  that  Cr-trapping  may  not  be  ’captured’  experimentally 
due  to  shortened  time  window  even  if  no  partitioning  of  Cr  occured  during 
nucleation  stage,  which  is  not  likely  to  occur  considering  that  the  lack  of 
partitioning  of  Cr  is  also  partially  attributed  to  its  lower  diffusivity. 

Enhanced  diffusivity  in  the  7/7'  interface  region  is  shown  to  operate  to  give  a 
substantial  increase  in  the  diffusion  kinetics  of  Al  over  what  would  be  expected 
from  mean  field  modeling.  While  the  concentration  profiles  in  the  interface 
region  are  largely  controlled  by  the  interplay  between  interfacial  energy  and 
chemical  bulk  energy,  the  bulk  composition  is  nevertheless  one  of  the  most 
important  factor  to  determine  the  actual  local  compositions  in  the  interface 
region.  Given  a  different  bulk  composition  as  we  studied  here,  the  degree  of 
enhancement  of  Al  vs.  Cr  diffuvisity  in  the  interface  region  can  be  decreased  or 
even  reversed.  An  example  of  such  is  shown  in  Fig.  9  for  the  mobility  profiles 
in  a  different  alloy,  i.e.  Ni-15Al-5Cr.  As  can  be  seen,  the  diffusivity  of  Al  in 
the  interface  region  is  significantly  decreased,  which  is  the  opposite  of  what  we 
observed  in  the  Ni-5.2Al-14.2Cr  alloy  as  shown  in  Fig.  5.  The  relative  mobility 
of  M^Cr/MAiAl  is  actually  increased  slightly  in  the  interface  region  as  shown 
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in  Fig.  9(c),(d).  These  results  indicate  that  the  Cr-trapping  phenomenon  is 
strongly  tied  to  the  selected  compositions  of  the  alloy. 

Phase-Field  modeling  indicates  that  the  Cr-trapping  can  only  be  observed 
for  non-equilibrium  critical  nucleus  concentrations.  Thermodynamic  modeling 
was  employed  to  show  that  non-equilibrium  critical  nucleus  concentrations 
were  possible.  Following  is  a  discussion  of  the  results  that  lead  us  to  making 
this  conclusion. 

T0  is  calculated  using  the  built-in  thermodynamic  database  and  the  result  is 
shown  in  Fig.  10.  In  the  chart,  a  shaded  zone  is  plotted  that  roughly  represents 
possible  nucleus  concentrations  that  can  survive.  The  triangle  zone  is  bounded 
by  the  T0  curve,  which  determines  if  the  nucleus  is  energetically  favorable,  the 
tie-line  and  the  line  with  no  partitioning  of  Cr.  The  tie-line  represents  an 
extreme  case  where  the  diffusivity  of  both  A1  and  Cr  are  so  high  that  it  is 
pure  thermodynamic  effect  in  controlling  nucleus  concentrational  profile.  On 
the  other  hand,  the  line  with  no  partitioning  of  Cr  represents  another  extreme 
case  where  the  diffusivity  of  Cr  is  so  small  that  it  is  totally  left  out  in  the 
formation  of  the  nucleus. 


4  Summary 


The  Cr-trapping  phenomenon  in  Ni-Al-Cr  alloys  is  studied  applying  a  recently 
developed  ternary  phase-field  model.  The  effect  of  the  compositional/ordering 
profiles  of  the  nucleus  and  the  diffusivity  on  the  evolution  of  Cr  concentration 
in  the  7'  precipitate  is  simulated.  The  simulation  work  reveals  that  the  Cr- 
trapping  phenomenon  can  only  be  reproduced  with  no  significant  partitioning 
of  Cr  at  the  nucleation  stage.  This  is  belived  to  be  a  major  factor  that  causes 
the  Cr-trapping  observed  by  experiments.  Another  factor  is  related  to  the 
much  stronger  composition-dependent  diffusivity  of  A1  as  compared  to  Cr  for 
the  alloy  composition  under  investigation.  The  diffuvisity  of  A1  is  enhanced 
by  approximately  an  order  of  magnitude  across  regions  that  have  significant 
concentration  inhomogeneity  at  a  7/7'  interface.  In  contrast,  the  diffusivity  of 
Cr  is  not  much  affected  in  the  interface  region.  This  is  another  major  factor 
that  has  caused  the  Cr-trapping  phenomenon. 
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Fig.  1.  Illustration  of  the  five  initial  structural/compositional  profiles  tested  in  this 
work. 
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(a) 


(b) 


(c)  (d) 

Fig.  2.  Effect  of  initial  nucleus  structural/compositional  profiles  on  Cr  concentration 
evolution  in  a  single  precipitate:  A  snapshot  of  the  precipitate  at  the  earlier  (a)  and 
later  stage  of  the  simulation  (b),  simulated  evolution  of  Cr  concentration  with  time 
(c)  and  expermentally  observed  Cr  concentration  evolution  (d). 
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Fig.  3.  Effect  of  initial  nucleus  structural/compositional  profiles  on  Cr  concentration 
evolution  in  multiple  precipitates. 


Fig.  4.  Mobility  effect  on  Cr  concentration  evolution  in  a  single  precipitate. 
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Fig.  5.  Profiles  of  the  mobility  (top  row)  and  the  ratio  of  M^Cr/MAllAl  (bottom 
row)  for  cases  (i)  and  (v)  along  the  middle  horizontal  cross-section:  (a)  and  (c) 
earlier  stage  at  1=0.2  (h);  (b)  and  (d)  later  stage  at  t=70  (h). 
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Fig.  6.  Profiles  of  the  gradient  of  the  thermodynamic  driving  force  for  cases  (i) 
and  (v)  shown  in  Fig.  2(b)  along  the  central  portion  of  the  middle  horizontal 
cross-section:  (a)  earlier  stage  at  t=0.2  (h),  (b)  later  stage  at  t=70  (h). 


distance  from  the  particle  center  distance  from  the  particle  center 


(a) 


(b) 


Fig.  7.  The  concentration  evolution  rate  profile  for  A1  and  Cr  for  cases  (i)  and  (v) 
shown  in  Fig.  2(b)  along  the  middle  horizontal  cross-section:  (a)  earlier  stage  at 
t=0.2  (h),  (b)  later  stage  at  t=70  (h). 
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Fig.  8.  Cr  concentration  evolution  in  a  single  precipitate  with  the  mobility  of  Cr 
set  to  equal  to  that  of  Al.  Reference  cases  are  included  for  comparison  where  the 
diffusion  mobilities  from  literature  database  are  used. 
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Fig.  9.  Profiles  of  the  mobility  (top  row)  and  the  ratio  of  M^Cr/MAllAl  (bottom 
row)  in  Ni-15Al-5Cr  alloys  along  the  middle  horizontal  cross-section:  (a)  and  (c) 
earlier  stage  at  t=0.2  (h);  (b)  and  (d)  later  stage  at  t=10  (h). 
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Fig.  10.  Calcualted  T0  curve  and  projected  survivable  nucleus  concentration  region. 
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